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ASYNCHRONOUS NETWORKS AND EVENT DRIVEN 
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Abstract. Real-world networks in technology, engineering and 
biology often exhibit dynamics that cannot be adequately repro¬ 
duced using network models given by smooth dynamical systems 
and a hxed network topology. Asynchronous networks give a the¬ 
oretical and conceptual framework for the study of network dy¬ 
namics where nodes can evolve independently of one another, be 
constrained, stop, and later restart, and where the interaction be¬ 
tween different components of the network may depend on time, 
state, and stochastic effects. This framework is sufficiently general 
to encompass a wide range of applications ranging from engineer¬ 
ing to neuroscience. Typically, dynamics is piecewise smooth and 
there are relationships with Filippov systems. In the first part of 
the paper, we give examples of asynchronous networks, and de¬ 
scribe the basic formalism and structure. In the second part, we 
make the notion of a functional asynchronous network rigorous, 
discuss the phenomenon of dynamical locks, and present a founda¬ 
tional result on the spatiotemporal factorization of the dynamics 
for a large class of functional asynchronous networks. 
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1. Introduction 

In this work we develop a theory of asynchronous networks and 
event driven dynamics. This theory constitutes an approach to net¬ 
work dynamics that takes account of features encountered in networks 
from modern technology, engineering, and biology, especially neuro¬ 
science. For these networks dynamics can involve a mix of distributed 
and decentralized control, adaptivity, event driven dynamics, switch¬ 
ing, varying network topology and hybrid dynamics (continuous and 
discrete). The associated network dynamics will generally only be 
piecewise smooth, nodes may stop and later restart and there may 
be no intrinsic global time (we give specihc examples and dehnitions 
later). Signihcantly, many of these networks have a function. For ex¬ 
ample, transportation networks bring people and goods from one point 






ASYNCHRONOUS NETWORKS 


3 


to another and neural networks may perform pattern recognition or 
computation. 

Given the success of network models based on smooth differential 
equations and methods based on statistical physics, thermodynamic 
formalism and averaging (which typically lead to smooth network dy¬ 
namics), it is not unreasonable to ask whether it is necessary to incor¬ 
porate issues such as nonsmoothness in a theory of network dynam¬ 
ics. While nonsmooth dynamics is more familiar in engineering than 
in physics, we argue below that ideas from engineering, control and 
nonsmooth dynamics apply to many classes of network and that nons¬ 
moothness often cannot be ignored in the analysis of network function. 
As part of these introductory comments, we also explain the motiva¬ 
tion underlying our work, and describe one of our main results: the 
modularization of dynamics theorem. 


Temporal averaging. Consider the analysis of a network where links 
are added and removed over time. Two extreme cases have been widely 
considered in the literature. If the network topology switches rapidly, 
relative to the time scale of the phenomenon being considered, then 
we may be able to replace the varying topology by the time-averaged 
topologji]. Providing that the network topology is not state dependent, 
the resulting dynamics will typically be smooth. On the other hand, if 
the topology changes slowly enough relative to the time scale of interest, 
we may regard the topology as constant and again we obtain smooth 
network dynamics. Either one of these approaches may be applicable 
in a system where time scales are clearly separated. 

However, in many situations, especially those involving control or 
close to bifurcation, changes in network topology may play a crucial 
role in network function and an averaging approach may fail or ne¬ 
glect essential structure. This is well-known for problems in optimized 
control where solutions are typically nonsmooth and averaging gives 
the wrong solutions (for example, in switching problems using ther¬ 
mostats). For an example with variable network topology, we cite the 
effects of changing connection structure (transmission line breakdown), 
or adding/subtracting a microgrid, on a power grid. Neither averag¬ 
ing nor the assumption of constant network structure are appropriate 
tools: we cannot average the problems away. Instead, we are forced to 
engage with an intermediate regime, where nonsmoothness (switching) 
and control play a crucial role in network function. 


^For example, if the input structure is additive - see section [5J 
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Spatial averaging and network evolution. Much current research 
on networks is related to the description and understanding of com¬ 
plex systems [TJ [181 EH El]. Roughly speaking, and avoiding a for¬ 
mal dehnition [H], we regard a complex system as a large network 
of nonlinearly interacting dynamical systems where there are feedback 
loops, multiple time and/or spatial scales, emergent behaviour, etc. 
One established approach to complex networks and systems uses ideas 
from statistical mechanics and thermodynamic formalism. For exam¬ 
ple, models of complex networks of interconnected neurons can some¬ 
times be described in terms of their information processing capability 
and entropy [SB]. These methods originate from applications to large 
interacting systems of particles in physics. As Schrodinger points out 
in his 1943 Trinity College, Dublin, lectures [59] 

“...the laws of physics and chemistry are statistical throughout.” 

In contrast to the laws of physics and chemistry, evolution plays a 
decisive role in the development of complex biological structure. Func¬ 
tional biological structures that provided the basis for evolutionary 
development can be quite small - the nematode worm caenorhabditis 
elegans has 302 neurons. If the underlying small-scale structure still 
has functional relevance, an approach based on statistical averages to 
complex biological networks has to be limited; on the one hand, averag¬ 
ing over the entire network will likely ignore any small scale structure, 
and on the other hand statistical averages have little or no meaning for 
small systems - at least on a short time scale. 

Reverse engineering large biological structures appears completely 
impractical; in part this is because of the role that evolution plays 
in the development of complex structure. Evolution works towards 
optimization of function, rather than simplicity, and is often local in 
character with the flavour of decentralized control. Similar issues arise 
in understanding evolved engineering structures. For example, the in¬ 
ternal combustion engine of a car in 1950 was a simple device, whose 
operation was synchronized through mechanical means. A modern in¬ 
ternal combustion engine is structurally complex and employs a mix 
of synchronous and asynchronous systems controlled by multiple com¬ 
puter processors, sensors and complex computer code. 

Reductionism. In nonlinear network dynamics, and complex systems 
generally, there is the question as to how far one can make use of re¬ 
ductionist techniques [5], [IH 2.5]. One approach, advanced by Alon 
and Kastan [39] in biology, has been the identihcation and descrip¬ 
tion of relatively simple and small dynamical units, such as non-linear 
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oscillators or network motifs (small network configurations that occur 
frequently in large biological networks [TSl Chapter 19]). Their premise 
is that a modular, or engineering, approach to network dynamics is fea¬ 
sible: identify building blocks, connect together to form networks and 
then describe dynamical properties of the resulting network in terms 
of the dynamics of its components. 

“Ideally, we would like to understand the dynamics of 
the entire network based on the dynamics of the indi¬ 
vidual building blocks.” Alon H page 27]. 

While such a reductionist approach works well in linear systems theory, 
where a superposition principle holds, or in the study of synchroniza¬ 
tion in weakly coupled systems of nonlinear approximately identical 
oscillators [55115SI El E2], it is usually unrealistic in the study of het¬ 
erogenous networks modelled by a system of analytic nonlinear differ¬ 
ential equations: network dynamics may bear little or no relationship 
to the intrinsic (uncoupled) dynamics of nodes. 

A theory of asynchronous networks. The theory of asynchronous 
networks we develop provides an approach to the analysis of dynamics 
and function in complex networks. We illustrate the setting for our 
main result with a simple example. Figure [T] shows the schematics of 
a network where there is only intermittent connection between node^. 
We assume eight nodes Ni,..., Ng. Each node W will be given an 
initial state and started at time Ti > 0. Crucially, we assume the 
network has a function: reaching designated terminal states in hnite 
time - indicated on the right hand side of the hgure. Nodes interact 
depending on their state. For example, referring to figure [H nodes W, 
W hrst interact during the event indicated by E“. Observe there is 
no global time dehned for this system but there is a partially ordered 
temporal structure: event always occurs after event E“ but may 
occur before or after event E^. We caution that while the direction of 
time is from left-to-right, there is no requirement of moving from left to 
right in the spatial variables: the phase space dimension for nodes could 
be greater than one and the initialization and terminations sets could 
be the same. This example can be generalized to allow for changes 
in the number and type of nodes after each event. The intermittent 
connection structure we use may be viewed as an extension of the idea 


^Figure [T] can be viewed as representing part of a threaded computer program. 
The events E“,..., will represent synchronization events - evolution of associ¬ 
ated threads is stopped until each thread has finished its computation and then 
variables are synchronized across the threads. 
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of conditional action as defined by Holland in the context of complex 
adaptive systems [3S]. 
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Figure 1. A functional feedforward network with 8 nodes 

Our main result, stated and proved in part II of this work |T3], is a 
modularization of dynamics theorem that yields a functional decompo¬ 
sition for a large class of asynchronous networks. Specifically, we give 
general conditions that enable us to represent a large class of functional 
asynchronous networks as feedforward functional networks of the type 
illustrated in hgure [H As a consequence, the function of the original 
network can be expressed explicitly in terms of uncoupled node dynam¬ 
ics and event function. Nonsmooth effects, such as changes in network 
topology through decoupling of nodes and stopping and restarting of 
nodes, are one of the crucial ingredients needed for this result. In net¬ 
works modelled by smooth dynamical systems, all nodes are effectively 
coupled to each other at all times and information propagates instantly 
across the entire network. Thus, a spatiotemporal decomposition is 
only possible if the network dynamics is nonsmooth and (subsets of) 
nodes are allowed to evolve independently of each other for periods of 
time. This allows the identihcation of dynamical units, each with its 
own function, that together comprise the dynamics and function of the 
entire network. The result highlights a drawback of averaging over a 
network: the loss of information about the individual functional units, 
and their temporal relations, that yield network function. 

A functional decomposition is natural from an evolutionary point of 
view: the goal of an evolutionary process is optimization of (network) 
function. Thus, rather than asking how network dynamics can be un¬ 
derstood in terms of the dynamics of constituent subnetworks - the 
classical reductionist question - the issue is how network function can 
be understood in terms of the function of network constituents. Our 
result not only gives a satisfactory answer to Alon’s question for a large 
class of functional asynchronous networks but suggests an approach to 
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determining key structural features of components of a complex sys¬ 
tem that is partly based on an evolutionary model for development of 
structure. Starting with a small well understood model, such as the 
class of functional feedforward networks described above, we propose 
looking at bifurcation in the context of optimising a network function 
~ for example, understanding the effect on function when we break the 
feedforward structure by adding feedback loops. 

Relations with distributed networks. An underlying theme and 
guide for our formulation and theory of asynchronous networks is that 
of efficiency and cost in large distributed networks. We recall the guide¬ 
lines given by Tannenbaum & van Steen [611 page 11] for scalability in 
large distributed networks (italicised comments added): 

• No machine has complete information about the (overall) system 
state, [communication limited] 

• Machines make decisions based only on local information, [decen¬ 
tralized control] 

• Failure of one machine does not ruin the algorithm, [redundancy] 

• There is no implicit assumption of global time. 

Of course, networks dynamics, in either technology, engineering or 
biology, is likely to involve a complex mix of synchronous and asyn¬ 
chronous components. In particular, timing (clocks, whether local or 
global) may be used to trigger the onset of events or processes as part 
of a weak mechanism for centralized control or resetting. Evolution 
is opportunistic - whatever works well will be adopted (and adapted) 
whether synchronous or asynchronous in character. In specihc cases, 
especially in biology, it may be a matter of debate as to which view¬ 
point - synchronous or asynchronous - is the most appropriate. The 
framework we develop is sufficiently flexible to allow for a wide mix of 
synchronous and asynchronous structure at the global or local level. 


Past work. Mathematically speaking, much of what we say has sig- 
nihcant overlap with other areas and past work. We cite in particular, 
the general area of nonsmooth dynamics, Filippov systems and hybrid 
systems (for example, [23 El EDI IH] ) and time dependent network struc¬ 
tures (for example, [3113 EH EZ])- While the theory of nonsmooth dy¬ 
namics focuses on problems in control, impact, and engineering, rather 
than networks, there is signihcant work studying bifurcation (for exam¬ 
ple [111 ED]) which is likely to apply to parts of the theory we describe. 
From a vast literature on networks and dynamics, we cite Newman’s 
text [52] for a comprehensive introduction to networks, and the very 
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recent tutorial of Porter & Gleeson IS7I which addresses questions re¬ 
lated to our work, gives an overview and introduction to dynamics on 
networks, and includes an extensive bibliography of past work. 

Outline of contents. After preliminaries in section [2l we give in sec¬ 
tion [3] vignettes (no technical details) of several asynchronous networks 
from technology, engineering, transport and neuroscience. In section HJ 
we give a mathematical formulation of an asynchronous network with a 
focus on event driven dynamics, and constraints. We follow in section 0 
with two more detailed examples of asynchronous networks including 
an illuminating and simple example of a transport network which re¬ 
quires minimal technical background yet exhibits many characteristic 
features of an asynchronous network, and a discussion of power grid 
models that indicates both the limitations and possibilities of our ap¬ 
proach. We conclude with a discussion of products of asynchronous 
networks in section E] that illuminates some of the subtle features of 
the event map. In part II [13], we develop the theory of functional 
asynchronous networks and give the statement and proof of the mod¬ 
ularization of dynamics theorem. 

Dedication. The genesis of this paper lies in a visit in 2010 by one us 
(MF) to work with Dave Broomhead at Manchester University. Dave 
was very interested in asynchronous processes and local clocks. During 
the visit, he came up with a 2 cell random dynamical systems model 
for the investigation of asynchronous dynamics and local time. This 
2 cell model provided the seed and stimulus for the work described in 
this paper. Dave’s illness and untimely death sadly meant that our 
planned collaboration on this work could not be realized. 

2. Preliminaries and generalities on networks 

2.1. Notational conventions. We recall a few mostly standard nota- 
tional conventions used throughout. Let N denote the natural numbers 
(the strictly positive integers), Z_|_ denote the set of nonnegative inte¬ 
gers, M+ = {x G M I X > 0}, and M(> 0) = {x G M+ | x 7^ 0}. Given 
n G N, dehne n = {1,..., n}. Let n* = {0,1, ..., n} and, more gener¬ 
ally, for A C N dehne A* = AU {0}. 

2.2. Network notation. We establish our conventions on network 
notation; we follow these throughout this work. 

Suppose the network Af has k nodes, W, • • •, Abusing notation, 
we often let Af denote both network and the set of nodes {A^i,..., N}f\. 


ASYNCHRONOUS NETWORKS 


9 


Denote the state or phase space for iVj by mJI and set M = Hiek ^ 
the network phase space. Denote the state of node iVj by x* G Mj and 
the network state by X = (xi,..., x^) G M. 

Smooth dynamics on M will be given by a system of ordinary differ¬ 
ential equations (ODEs) of the form 

(1) x'=/i(xi;xj,,...,xj-^,), iGk, 

where the components ft are at least (usually C°° or analytic) and 
the following conditions are satished. 

(Nl) For alH G k, ji < ... < are distinct elements of k \ {i} (and 
so Ci < k). 

Set J(i) = {ji,..., jej C k, i G k (J(i) may be empty). 

(N2) For each z G k, the evolution of iVj depends nontrivially on the 
state of Nj, j G J{i), in the sense that there exists a choice of x* G Mi 
and Xj^ G jg G ^(f) \ {j}, such that /j(xj; x^j,..., Xj^J is not 
constant as a function of Xj. 

(N3) We generally assume the evolution of W depends on the state of 
Ni. If we need to emphasize that fi does not depend on x* in the sense 
of (N2), we write /i(xj^,... if J{i) 7 ^ 0- If J{i) = 0, we regard 
the dependence of fi on x* as nontrivial iff fi is not identically zero and 
then write /i(xj). Otherwise /j = 0. 

Remark 2.1. Given network equations ([T]) which do not satisfy (Nl- 
3), we can Erst redehne the fi so as to satisfy (Nl). Next we remove 
trivial dependencies so as to satisfy (N2). Finally, we check for the 
dependence of fi on the internal state Xj and modify the fi as necessary 
to achieve (N3). If fi = 0, we can remove the node from the network. 
Consequently, it is no loss of generality to always assume that (Nl-3) 
are satished, with /j ^ 0. A consequence is that any network vector 
field f = (/i,...,/fc) : M —TM can be uniquely written in the 
form ([T]) so as to satisfy (Nl-3). 

Let M{k) denote the space of k x k 0-1 matrices (d = (/3ij)jjgk with 
coefficients in {0,1} and (da = 0, all z G k. Each fd G M{k) determines 
uniquely a directed graph with vertices Ni,...,Nk and directed 
edge Nj —)■ W iff /dij = 1 and i ^ j. The matrix (d is the adjacency 
matrix of T We refer to /9 as a connection structure on M. 

^We assume the phase space for each node is a connected differential manifold - 
usually a domain in R" or the n-torus, T". 
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If f : M —)■ TM is a network vector field satisfying (Nl-3), then f 
determines a unique connection structure C(f) G M{k) with associated 
graph r( 7 (f) • In order to specify the graph uniquely, it suffices to specify 
the set of directed edges. 

We dehne the network graph T = r(A/', f) to be the directed graph 
r, 7 (f). Thus, r(7\/', f) has node set Af = {Ni,... ,Nk} and a directed 
connection Nj —)• W will be an edge of T if and only if j ^ i and the 
dynamical evolution of W depends nontrivially on the state of Nj. 

Remark 2.2. Our conventions are different from formalisms involving 
multiple edge types (for example, see [321 E] for continuous dynamics 
and [1] for discrete dynamics). We allow at most one connection be¬ 
tween distinct nodes of the network graph and do not use self-loops: 
connections encode dependence. 

2.2.1. Additive input structure. In many cases of interest, we have an 
additive input structure [26] and the components fi of f may be written 

Si 

(2) /i(xi; ... , XjJ = Fi(xi) -f %(xj^, Xi), i G k. 

S=1 

Additive input structure implies that there are no interactions between 
inputs Nj, Nk Ni, as long as j, k ^ i, j ^ k, and allows us to add 
and subtract inputs and nodes in a consistent way. We may think of 
x' = Fj(xj) as dehning the intrinsic dynamics of the node. 

Remarks 2.3. (1) Additive input structure is usually assumed for mod¬ 
elling weakly coupled nonlinear oscillators and is required for reduction 
to the standard Kuramoto phase oscillator model [l2l [2^ [38] . 

(2) If we identify a null state z* for each node Nj, then the decompo¬ 
sition ([2D will be unique if we require Fjj(z*,Xj) = C@. If a node is in 
the null state then it has no output to other nodes and is ‘invisible’ to 
the rest of the network. If we have an additive structure on the phase 
spaces Mi (for example, each Mi is a domain in or an n-torus T") 
it is natural to take z* = 0. 

(3) If Mi = or T”, i G k, and Fjj^(xj^,Xj) = Gjj,(xj^ -x*), i,js G k, 
the coupling is diffusive (see [U §2.5] for general phase spaces). 

2.3. Synchronous networks. Systems of ordinary differential equa¬ 
tions like ([ID give mathematical models for synchronous networks. By 
synchronous, we mean nodes are all synchronized to a global clock - the 
terminology comes from computer science. Indeed, if each node comes 

^For identical phase spaces, assume inputs are asymmetric - Fij Fn, if j ^ t. 
For symmetric inputs see [SHE]. 
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with a local clock, then all the clocks can be set to the same time pro¬ 
vided that the network is connected (we ignore the issue of delays, but 
see Ha). The synchronization of local clocks is essentially forced by the 
model and the connectivity of the network graph; nodes cannot evolve 
independently of one another unless the network is disconnected. 

We recall some characteristic features of synchronous networks. 

Global evolution: Nodes never evolve independently of each 
other: if the state of any node is perturbed, then generically 
the evolution of the states of the remaining nodes changes. 

Stopped nodes: If a node (or subset of node variables) is at equi¬ 
librium or “stopped” for a period of time, it will remain stopped 
for all future time. If a node has a non-zero initialization, it will 
never stop (in hnite time). 

Fixed connection strnctnre: The connection structure of a syn¬ 
chronous network is hxed: it does not vary in time and is not 
dependent on node states: one system of ODEs suffices to model 
network dynamics. 

Reversibility: Solutions are uniquely dehned in backward time. 

3. Asynchronous networks: examples 

In this section, we give several vignettes of asynchronous networks 
that illustrate the main features differentiating them from synchronous 
networks. We amplify two of these examples in section [5] after we have 
developed our basic formalism for asynchronous networks. 

Example 3.1 (Threaded and parallel computation). Threaded or par¬ 
allelized computation provides an example of a discrete stochastic asyn¬ 
chronous network. Computation based on a single processor (or sin¬ 
gle core of a processor) proceeds synchronously and sequentially. The 
speed of the computation is dependent on the clock speed of the pro¬ 
cessor as the processor clock synchronizes the various steps in the com¬ 
putation. In threaded or parallel computation, computation is broken 
into blocks or ‘threads’ which are then computed independently of each 
other at a rate that is partly dependent on the clock rates of the proces¬ 
sors involved in the computation (these need not be identical). At cer¬ 
tain points in the computation, threads need to exchange information 
with other threads. This process involves stopping and synchronizing 
(updating) the thread states: a thread may have to stop and wait for 
other threads to complete their computations and update data before 
it can continue with its own computation. 

Threaded computation is non-deterministic: the running and stop¬ 
ping times of each thread are unpredictable and differ from run to run. 
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Each thread has its own clock (determined by its associated proces¬ 
sor). Threads will be nnaware of the clock times of other threads except 
dnring the stopping and synchronization events which can be managed 
synchronously (central control) or asynchronously (local control). 

This example shows many characteristic features of an asynchronous 
network: nodes (threads) evolving independently of each other, and 
stopping, synchronization and restarting events. The network also has 
a function - transforming a set of initial data into a set of hnal data 
in hnite time - and there is the possibility of incorrect code that can 
lead to a process that stops before the computation is complete (a 
deadlock), or errors where threads try to access a resource at the same 
time (race condition). <0 

Example 3.2 (Power grids & microgrids). A power grid consists of a 
connected network of various types of generators and loads connected 
by transmission lines. A critical issue for the stability of the power 
grid is maintaining tight voltage frequency synchronization across the 
grid in the presence of voltage phase differences between generators and 
loads and variation in generator outputs and loads. We refer to Kun- 
dur 110] for classical power grid theory, Dorfler et al. [23] or Nishikawa 
& Motter [53], for some more recent and mathematical perspectives, 
and [H] for general issues and dehnitions on power system stability. 

Historically, power grids have been centrally controlled and one of 
the main stability issues has been the effect on stability of a sudden 
change in structure - such as the removal of a transmission line, break¬ 
down of a generator or big change in load. Detailed models of the 
power grid need to account for a complex multi-timescale stiff system. 
Typically stability has been analyzed using numerical methods. How¬ 
ever, relatively simple classes of network models for power grids based 
on frequency and phase synchronization have been developed which are 
applicable for the analysis of some stability and control issues, espe¬ 
cially those described in the next paragraph. We describe these models 
in more technical detail in section [5l 

Interest has recently focused on renewable (small) energy sources in 
a power grid (for example, wind and solar power) and how to inte¬ 
grate microgrids based on renewable sources into the power grid using 
a mix of centralized and decentralized control. Concurrent with this 
interest is the issue of smart grids: modifying local loads in terms of 
the availability and real time costs of power. While the classical power 
grid model is of a synchronous network, though with asynchronous fea¬ 
tures such as the effects on stability of the breakdown of a connection 
(transmission line), these problems focus on asynchronous networks. 
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For example, given a microgrid with renewable energy sources such as 
wind and solar, time varying loads and buffers (large capacity batter¬ 
ies), how can the microgrid be switched in and out of the main power 
grid while maintaining overall system stability? In this case, switch¬ 
ing will be determined by state (for example, frequency changes in the 
main power grid signifying changes in power demand or changes in the 
output of renewable sources or battery reserves) and stochastic effects 
(resulting, for example, from load changes and the incorporation of 
smart grid technology). This is already a tricky problem of distributed 
and decentralized control with just one microgrid; in the presence of 
many microgrids there is the potential problem of synchronization of 
switching microgrids in and out of the main power grid. Similar prob¬ 
lems occur in smart grids |62j . 

Asynchronous features of power grid networks include variation in 
connection and node structure (separation, or islanding, of microgrids 
from main power grid), state dependence of connection structure, syn¬ 
chronization and restarting events (during incorporation of microgrid 
into main grid). <0 

Example 3.3 (Thresholds, spiking networks and adaptation). Many 
mathematical models from engineering and biology incorporate thresh¬ 
olds. For networks, when a node attains a threshold, there are often 
changes (addition, deletion, weights) in connections to another nodes. 
For networks of neurons, reaching a threshold can result in a neuron 
firing (spiking) and short term connections to other neurons (for trans¬ 
mission of the spike). For learning mechanisms, such as Spike-Timing 
Dependent Plasticity (STDP) [2H] relative timings (the order of firing) 
are crucial [301 [m EH and so each neuron, or connection between a 
pair of neurons, comes with a ‘local clock’ that governs the adaptation 
in STDP. In general, networks with thresholds, spiking and adapta¬ 
tion provide characteristic examples of asynchronous networks where 
dynamics is piecewise smooth and hybrid - a mix of continuous and 
discrete dynamics. Spiking networks also highlight the importance of 
efficient communication in large networks: spiking induced connections 
between neurons are brief and low cost. There is also no oscillator clock 
governing all computations along the lines of a single processor com¬ 
puter. These examples all fit well into the framework of asynchronous 
networks but, on account of the background knowledge required, we 
develop the theory and formalism elsewhere [T3|. <) 

Example 3.4 (Transport & production networks). We discuss trans¬ 
port networks first. For simplicity, we work with a single transport 
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mode: trains. Typically, trains have to be schednled to be in a sta¬ 
tion for overlapping times (stopping, restarting, connections and local 
times) so that passengers can transfer between trains, or stop in a pass¬ 
ing loop (so that trains can pass on a single track line). In addition, a 
train can divide into two parts or two trains can be combined (varia¬ 
tion in node strnctnre, stopping and synchronization event). Generally, 
transport networks will have asynchronons featnres and exhibit state 
dependent connection structure, local times and have a strong stochas¬ 
tic component (for example, in stopping and restarting times). We de¬ 
velop a simple formal transport model in section O (and in part II [T3] ) 
that illustrates basic ideas and results in the theory of asynchronous 
networks but does not require extensive background knowledge. 

A simple example of a production network is a paint mixer. Assume 
a controller which accepts inputs - requested colour - which, after 
computation to hnd tint weights (‘tint code’), signals a request to inject 
selected tints according to the tint code into the base paint which is 
then mixed. The output is a can of coloured and fully mixed paint. 
Dynamics plays a limited role - except possibly at the mixing stage 
(for example, if there is a sensor that can detect an acceptable level 
of mixing). For this network, there is a varying connection structure 
determined by the signalling and tint injection. A characteristic feature 
of this, and many production networks, is the large variation in time 
scales. Signalling will typically be very fast, injection moderately fast 
and mixing rather slow. If the times of inputs to the controller are 
stochastic (for example, follow a Poisson process), then there will be 
issues of queueing and prioritization of inputs. If it is intended to 
maximize usage of the production facilities and avoid long waits, then it 
is natural to suppose that there are several mixing units and the output 
of the tint units is switched between mixer units according to their 
availability. Of course, the paint mixer may be a small part of a much 
larger distributed production network for which we can expect multiple 
time scales, switching between production units, changing the output 
of production units, stopping or restarting units, etc. The control of 
large distributed production systems will typically involve a mix of 
decentralized and centralized control. 

Synthesis of proteins at the cellular level can be viewed as a gener¬ 
alization of the paint mixer model. We refer the reader to Alon [U 8, 
Chapter 1] for background and more details, especially on transcription 
networks. <0 

We summarize some of the key features of asynchronous networks 
illustrated by all of the preceeding examples. 
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(1) Variable connection structure and dependencies between nodes. 
Changes in connection structure may depend on the state of the 
system or be given by a stochastic process. 

(2) Synchronization events associated with stopping or waiting states 
of nodes. 

(3) Order of events may depend on the initialization of the system 
or stochastic effects. 

(4) Dynamics is only piecewise smooth and there may be a mix of 
continuous and discrete dynamics. 

(5) Aspects involving function, adaptation and control. 

(6) Evolution only defined for forward time - systems are non- 
reversible. 

4. A Mathematical model for asynchronous networks 

In this section we formalize the notion of an asynchronous network. 
Our focus is on deterministic (not stochastic) and continuous time asyn¬ 
chronous networks which are autonomous (no explicit dependencies on 
time) and we use the term ‘asynchronous network’ as synonym for a 
deterministic and autonomous continuous time asynchronous network. 

4.1. Basic formalism for asynchronous networks. Consider a net¬ 
work M with k nodes, Ni,..., Nk, and follow the conventions of sec¬ 
tion |2j each node Ni has phase space Mj, and M = ~ 

network phase space. A network vector field f on M is assumed to sat¬ 
isfy conditions (Nl-3) and so determines a unique connection structure 
C(f) G M{k) and associated network graph rc'(f) (no self-loops). 

Stopping, waiting, and synchronization are characteristic features of 
asynchronous networks. If nodes of a network are stopped or partially 
stopped, then node dynamics will be constrained to subsets of node 
phase space. We codify this situation by introducing a constraining 
node Nq that, when connected to W, implies that dynamics on W is 
constrained. We give the precise dehnition of constraint shortly (in 
4.3); for the present, the reader may regard a constrained node as 
stopped - node dynamics is defined by the zero vector field. We only 
allow connections Nq ^ Ni, i E ]<., and do not consider connections 
Ni — )■ Nq, i G k*. Henceforth we usually always assume there is a 
constraining node and let Af = {Nq, Ni,..., N^} denote the set of 
nodes. We emphasize that the constraining node Nq has no dynamics 
and no associated phase space. In a network with no constraints (there 
are no connections Nq —>■ W), the constraining node Nq plays no role 
and can be omitted. If we allow constraints, there may be more than 
one type of constraint on a node A,. 
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Suppose that there are Pi E N different constraints on the node iVj, 
i E k. Set P = (pi,... ,pk) G and let M,{k] P) denote the space of 

k X {k + 1) matrices a = (cipOiekjGk* such that 

(1) e M{k) (and so an = 0,i G k). 

(2) ciio G p*, f G k. 

If q; G M,{k] P), we dehne the directed graph Fq, by 

(1 ) Fq, has node set A/”. 

(2) For all i, j G k, Nj —)■ iVj is an edge iff = 1. 

(3) Nq -e Ni is an edge iff aio 7 ^ 0. We write Nq -a iV* if we need 

to specify the constraint corresponding to £ G p*. 

We usually abbreviate M,(/c;P) to M,{k). Let 0 G M,{k) denote the 
empty connection structure (no edges). 

If a G M,{k), let a^ denote the hrst column (Q;io)jGk of a. We have 
a natural projection tt : M,{k) -E M{k); a 1 —)■ a*', defined by omitting 
the column a^. We write a E M,{k) uniquely as 

a = (q;° I a*'). 

The column vector codihes the connections from the constraining 
node and ct encodes the connections between the nodes {Ni ,..., N^}. 

Let a E M,{k). We provisionally dehne an a-admissible vector field 
f = (/i,..., /fc) to be a network vector held such that for i,j G k, i 7 ^ j, 
fi depends on the state of Nj ih = 1. If there is a connection 
Nq -e Ni [aio 7 ^ 0), then there is a nontrivial constraint on W- An 
a-admissible vector held has constrained dynamics if there are connec¬ 
tions from the constraining node. If a = 0, nodes are uncoupled and 
unconstrained. 

Definition 4.1. (Notation and assumptions as above.) 

(1) A generalized connection structure Al is a (nonempty) set of 
connection structures on Af. 

(2) An A-structure is a set = {f" | a G A} of network vector 
helds such that each f" G is a-admissible. 

Interactions between nodes in asynchronous networks may vary and 
can be state or time dependent or both. We focus on state dependence 
and assume interactions and constraints are determined by the state of 
the network through an event map T : M —)■ Al. 

Definition 4.2. Given a network Af, generalized connection structure 
A, Al-structure iF, and surjective event map T ; M —)■ Al, the quadruple 
91 = {Af,A,iF,S) dehnes an asynchronous network 
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The network vector field of T1 is given by the state dependent vector 
field F ; M ^ TM defined by 

F(X) = X G M. 

Remarks 4.3. (1) Subject to simple regularity conditions, which we give 
later, the network vector field F will have a uniquely defined semifiow. 
(2) In the sequel we often use the notation ffi as shorthand for the 
asynchronous network (A/*, A,J^,S) (by extension, 01“ will be shorthand 
for {N‘°‘, A°‘, , 8 °'), etc.). 

Example 4.4. Let k = 2 and Mi = M 2 = M x T. Suppose that 
dynamics of the uncoupled node iVj is given by the smooth vector field 
Vi{xi, 6 i) = where fi{ 0 ) 7 ^ 0 , ca* G M, i G 2 . 

Assume constrained dynamics for either node is defined on the invari¬ 
ant circle {0} x T C M x T by the vector field Zi{xi, 9i) = (0, a;*), i G 2. 
When both nodes are constrained {xi = 0:2 = 0), assume (constrained) 
coupling is defined by the vector field H = (iLi, H 2 ), where 

//i(a;i, 6^1, a;2, ^ 2 ) = (0, Wi + h(02 ~ ^ 1 )) 

7/2(2^!, ^1, a; 2 , ^2) = ( 0 , Ci ;2 + h( 0 i — ^2)), 

and h : T —)■ M is smooth. The 2-tori {{xi,X 2 )} x are invariant by 
the flow of H for all {xi,X 2 ) G Revert to standard (uncoupled and 
unconstrained) dynamics when |^i — ^ 2 ! < where 0 < e 1. We 
describe the network dynamics using asynchronous network formalism. 

Take the generalized connection structure A = {0, ai, 02, /9}, where 
Q.i = TVq —t fVj, i G 2, and (3 = Nq —v Ni GG N 2 i — Nq. 

Take = {H | 7 G A}, where 

f® = (I/i,F 2), = (^1,R2), f“^ = (Ri,^ 2), f"" = {Hi,H2). 

Define the event map : (M x T)^ —)■ Al by 

T( 0 , 6 ^ 1 , 0 , 62 ) = (3, if \di — 021 > £ 

= 0, if \9i — 021 < £ 

£{0,9i,X2,92) = tti, if a ;2 7 ^ 0 
S{xi,9i,0,92) = 0 : 2 , if 7 ^ 0 
S(xi,9i,X2,92) = 0 , if a;ia ;2 7 ^ 0 . 

Network dynamics is given by the vector field F(X) = f^*^^)(X). Tra¬ 
jectories for F are built from pieces of the trajectories of f®, f"b 
and f^. Using the condition /i(0) 7 ^ 0, i G 2, we see easily that F has 
a well-defined semifiow <ht(a;i, 6 ^ 1 , 0 : 2 , 6 * 2 ), which is continuous in time 
t > 0 but is not necessarily continuous in (xi, 6 * 1 , 0 : 2 , 6 * 2 ). <0 
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4.2. Local foliations. Conditions for a constrained node iVj will be 
given in terms of foliations of open subsets of Mj. We start by recalling 
basic definitions on foliations (see [16] for a detailed review). 

A p-dimensional smooth (always C°° here) foliation C of the m- 
dimensional manifold W consists of a partition {L^ | a G A} of hh 
into connected sets, called leaves, such that for every x G W, we can 
choose an open neighbourhood U oi x and smooth embedding ip : 
U such that for each leaf La, the components of (j){La fl U) are 

given by equations x^^^ = constant,..., x"^ = constant. Each leaf of a 
foliation will be an immersed p-dimensional submanifold of W. For our 
applications, we always assume leaves are properly embedded closed 
submanifolds of W, p < m, and that the manifold W has finitely many 
connected components. In general, a smooth foliation of the manifold 
W will consist of a smooth foliation of each connected component of 
W such that the dimension of leaves is constant on each connected 
component of W. 

Examples 4.5. (1) Every smooth nonsingular vector held on W de- 
hnes a 1-dimensional smooth foliation of W (“how-box” theorem of 
dynamical systems). The leaves are trajectories of the vector held. 

(2) If hh = A X B, where A and B are manifolds, we have the product 
foliations C{A) and C{B) of W dehned by C{A) = {A x {6} | 6 G S} 
and C{B) = {{a} x i? | a G A}. Each leaf C{A) is transverse to every 
leaf of C{B). More generally, foliations C, C are transverse if leaves are 
transverse. A foliation of W, even by compact 1-dimensional leaves, 
need not have a transverse foliation. The best-known example is the 
Hopf hbration which dehnes a foliation of S'^ into circles. <0 

Suppose that £ is a p-dimensional smooth foliation of W with leaves 
{La I a G A}. The tangent bundle along the foliation r : L —)■ hh is the 
smooth vector sub-bundle of the tangent bundle TW of W dehned by 

L = IJ C TW. 

X£La , CtGA 

4.3. Constrained nodes and admissible vector fields. Following 

section im we assume Af = [Nq, Ni,..., N^}, where the nodes W have 
phase space Mj, i G k. Fix a fc-tuple P = (pi,... ,pk) G In what 
follows, we assume P 7 ^ 0. 

Definition 4.6. (Notation and assumptions as above.) A family C = 
{(Wj, Ci) I i G k} is a constraint structure on M if, for alH G k with 
Pi > 0, 

(1) Wj = {Wf I £ G Pi} is a family of nonempty open subsets of 
Mj. 
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(2) Ci = 11 e Pi}, where C\ is a smooth foliation of hh/. 

Remarks 4.7. (1) If Pi = 0, there are no constraints on iVj. 

(2) If Pi = 1, we set Wj = (14}, £j) and Ci is a smooth foliation of the 
nonempty open subset 14} of Mj. If we allow the dimension of leaves 
to vary between different connected components, and the families Wj 
to consist of disjoint open subsets of Mj, z e k, then we can reduce 
to the case Pj < 1 by taking 14} = IJ^ W/ and £j to be the foliation 
determined on 14} by Ci\Wf = £}, £ G pi. For our applications, it is no 
loss of generality to assume that Wj always consists of disjoint open 
subsets of Mj, z e k. 

We can now give a precise definition of an a-admissible vector field 
when there are constraints. 

Definition 4.8. Fix a constraint structure C = {(Wj,/7j) | z G k} on 
M and let a G M,{k). A smooth vector field f = (/i,..., on M is 
an a-admissible vector field if 

(1) For z, j G k, z 7 ^ j, fi depends on iff Oj^ = 1. 

(2) If Ojo = £ > 0, then /j is tangent to the smooth foliation 

at all points of IF/ C Mj. Equivalently, ffWf defines a section 
of L}, the tangent bundle along the foliation Cf. 

Example 4.9. Suppose that pj = 1 and Ojo = 1 so that there is a 
constraining connection Nq —>■ A}- Let f = (/i,..., /a,) be a-admissible. 
Mi = R^, and Ci be an (£ — p)-dimensional foliation of Mj with leaves 
given hy Xr^ = ci,..., Xr^ = Cp. The components //S ..., f-"’ of /j = 
(//,..., //) will be identically zero and the node A} is partially stopped 
on each leaf. This is the situation described in example 14.41 where the 
1-dimensional foliation of R x T is {{x} x T | a; G R}. <0 

Remark 4.10. Note that if A"o —)■ A} •(— A}-, then the coupling from A}- 
must respect constraints on A} though now of course the dynamics on 
a leaf of Ci will depend on the state of A}. 

4.4. The event map. Let .4. be a generalized connection structure 
with constraint structure C = {(Wj, Ci) \ i G k}. Let : M —)■ 4, be 
an event map and recall S is always assumed to be surjective. 

For each a G 4., define the event set E" C M by 

= {X G M|T(X) = a}. 

The event sets (E" | a G A} partition the network phase space M. 
We require additional conditions on the event map when there are 
constraints. These conditions relate the event sets to the constraint 
structure C and are required because foliations are only locally defined. 
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Let TTj : M —)• Mi denote the projection map onto the phase space of 
iVj, z G k. Given z G k, £ G Pi, dehne 

El= IJ 7r,(E“) C M,. 

{a I aio=e} 

Definition 4.11. The event map : M —)■ .4, is constraint regular if 
for all z G k, £ G Pi, we have 

Henceforth we assume that event maps are constraint regular. 

4.5. Asynchronons network with constraints. 

Definition 4.12. An asynchronous network Tl = {J\f , A, E, S), with 
constraint structure C, consists of 

(1) A hnite set Af = {Aq, Ni, ..., N^} nodes with associated phase 
spaces Mj, z G k. 

(2) A generalized connection structure A C M,(k). 

(3) An A-structure T = {f" | a G A} consisting of admissible vec¬ 
tor fields. 

(4) A (constraint regular) event map £ \ A/l ^ A. 

Remark 4.13. If A consists of a single connection structure a (with 
or without constraints), then E consists of one vector held f = f“, 
with dependencies given by a. We recover a synchronous network with 
dynamics dehned by f and a hxed connection structure. 

4.6. Network vector field of an asynchronous network. An asyn¬ 
chronous network Tl uniquely determines the network vector field F by 

(3) F(X) = f^W(X), X G M. 

Remarks 4.14. (1) We may give a discrete version of dehnition 14.121 
each f“ will be a network map f" : M —)■ M and dynamics is dehned 
by the map F : M —)■ M given by ([3]). 

(2) Equation (|3]) dehnes a state dependent dynamical system. Similar 
structures have been used in engineering applications (for example, 
[54]). We indicate in section 15. 1.31 a relationship with Filippov systems 
(this is explored further in [12] )• However, the notion of an integral 
curve for an asynchronous network is generally diherent from that of a 
Filippov system, see examples I4.17f 2j. 

(3) The network vector held does not uniquely determine A, £ or E. 
Usually, however, the choice of A, £ and T is naturally determined by 
the problem. Sometimes it is convenient to view the network vector 
held as the basic object and regard asynchronous networks as being 
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equivalent if they define the same network vector field. 

(4) Since the event sets {E^ \ a G A\ partition M, the network vector 

field F only depends on Rather than assnme that f" is smooth 

on M, we could have required that each f“ was dehned as smooth map 
in the sense of Whitney [63] on (and so extends smoothly to M). 

(5) Although the vector fields f" G are assumed to satisfy (Nl-3), 

this may not hold for a E A. Sometimes, but not always, 

there is an equivalent network W such that the dependencies of each 
admissible vector held for W are not changed by restriction to the 
corresponding event set. 

4.7. Integral curves and proper asynchronous networks. We 

start with a dehnition of integral curve suitable for asynchronous net¬ 
works. 

Definition 4.15. Let fh be an asynchronous network with network 
vector held F. An integral curve or trajectory for F with initial condi¬ 
tion Xq G M is a map 4> : [0, T) —>■ M, T G (0, oo], satisfying 

(1) 0(0) = Xo. 

(2) (p is continuous. 

(3) There exists a closed countable subset D of [0,T) such that for 
every u E D, there exists v E D U {T}, v > u, such that 

(a) {u, v) r\ D = %. 

(b) 0 is on {u,v) and (p'if) = F(0(t)), t G {u,v). 

(c) limi^„+0'(t) =F(0 (m)). 

Remarks 4.16. (1) It is routine to verify that if 0 : [0, S') —?■ M 
is another integral curve with initial condition Xq, then 0 = 0 on 
[0, minis', T}) (uniqueness). As a consequence we can dehne the maxi¬ 
mal integral curve 0 : [0,Tinax) M with initial condition Xq. In the 
sequel, integral curves will be maximal unless otherwise indicated. 

(2) If T = oo in the dehnition, the trajectory 0 : M_|_ —)■ M is complete. 

(3) The set D may have accumulation points in D - accumulation is 
always from the left on account of condition (3a). In the examples we 
consider D will always be a hnite set. 

(4) Typically, for each u E D, there exists a E A such that S(0(t)) = a 
for t G {u,v) and so <p{{u,v)) C E°‘. Condition (3c) implies that if 
S{4>{u)) = /3 0 a, we must have f"(0(n)) = f^(0(n)). 

Without further conditions on the event map, the vector held F de¬ 
termined by an asynchronous network D1 may not have integral curves 
through every point of the phase space. 

Examples 4.17. (1) Take event sets E^ = {{xi,X 2 ) | xi < 0}, = 

\ E^, and corresponding constant vector helds = (1, —2), = 
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Figure 2. Integral curves for the network vector field 
may not be well defined (a) and may differ from those 
given by the Filippov conventions (b). 


(—1,0) (see hgure [2](a)). Trajectories cannot be continued, according 
to dehnition 14.151 once they meet xi = 0. One way round this problem 
is to dehne a new event set = dE^ and the sliding vector held 
+ f^ = (0, —2). There is then a complete integral curve through 
every point of and the corresponding semihow $ : x R+ —)■ is 

continuous. This approach is based on the Filippov construction [271 
Chapter 2, page 50] where we take a vector held in the positive cone 
dehned by f^, (often the unique convex combination Af^ + (1 ~ 
which is tangent to dE^ = E^). 

(2) Take event sets E^ = {{xi,X 2 ) \ Xi ^ X 2 }, E"^ = {(ti, X 2 ) \ Xi = X 2 }, 
and corresponding vector helds ^^{x\,X 2 ) = (1, —1), f^(xi,X 2 ) = (0,0) 
(see hgure [2](b) and note that the event E"^ models a collision, after 
which dynamics stops). Integral curves are dehned for all initial condi¬ 
tions in but the semihow $ : x R+ R^ will not be continuous 

on E^. Here the Filippov construction gives the wrong network solution 
- the diagonal is regarded as a removable singularity. 

We discuss the relationship between asynchronous networks and Fil¬ 
ippov systems further in section 15.1.31 see also [12] . <0 

Definition 4.18. The asynchronous network Tl is proper if for all 
X G M, the maximal integral curve through X is complete: ■ 

[0, cx)) —M. 

Remarks 4.19. (1) If Tl is proper, network dynamics is given by a 
semihow 4> : M x R+ —M. Although <I)(X, f) will be continuous as a 
function of f G R+, it need not be continuous as a function of X G M 
(see examples 14 .1 7f 21). 

(2) In many cases of interest, some of the node phase spaces Mj may be 
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open domains in M"" with with dMi ^ 0. Here there is the possibility 
that trajectories may exit M: if </> = , cf)^) is a trajectory, there 

may exist i G k and a smallest s > 0 such that (t>i{s) limi_j. 5 _ 0j(t) G 
dMi. The maximal domain for 0 is necessarily [0, s). Under additional 
hypotheses, it may be possible to extend 0 to a complete trajectory by 
setting Fj = 0 on M"' \ M^, j G k (the jth component of 4> is stopped 
when it meets the boundary of Mj). In this way, we can regard Dd as 
proper. We develop this point of view further in part II [T3] . 

Event sets are typically dehned by analytic and algebraic conditions 
that reflect logical conditions on the underlying dynamics. 

Definition 4.20. Let Dd be an asynchronous network. The event struc¬ 
ture {E® I a G .4.} of Dd is regular if the event sets are all semiana- 
lytic subset^ of M. 

Remark 4.21. For the examples in this paper, event sets will typically 
be semialgebraic - dehned by polynomial equalities and inequalities. 

Definition 4.22. An asynchronous network Dd is amenable if 

(1) The event structure {U" | a G A} is regular. 

(2) If X G U", a G 4., there exists a maximal t(X) G (0, cxo] such 
that the integral curve 0x through X is dehned on [0,t(X)) 
and 

tG[0,t(X)). 

(3) Either M* is compact without boundary or Mj = and vector 
helds have at most linear growth on Mj: 3a, 6 > 0 such that 

||f“(X)|| < a + 6||X|1, X G M, aeA. 

Remarks 4.23. (1) Condition (2) of dehnition 14.221 suggests that the 
vector held f" should in some sense be tangent to E". The issue of 
tangency can be made precise using the regularity assumption which 
implies that E" has a locally hnite stratihcation into submanifolds 
without boundary (for example, the canonical Whitney regular strat¬ 
ihcation of each event set imiiH]). This allows us to unambiguously 
dehne tangency at points of which do not lie in the boundary of 
strata. Care is needed at points lying in the boundary of strata and 
in the example below we indicate how the geometric structure of the 
event set can impose strong constraints on associated vector helds. 

(2) If an event set is a closed submanifold without boundary, it follows 
from dehnition I4.22f 2l that any trajectory that meets the event set will 

^Defined locally by analytic equations and inequalities. We refer to (JU [TS] for 
precise definitions and properties. 
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never leave the event set. 

(3) In part II we extend definition 14.221^ 3) to allow for trajectories to 
exit the domain and stop (see remark I4.19f 2lh 

(4) We may extend the dehnition of amenability to include asynchro¬ 
nous networks which are equivalent to an amenable network. 

Examples 4.24. Take k = 2, Mi = M 2 = M. 

(1) As event sets take the semialgebraic subsets of dehned by 

2 

Ei = {(a;, 0 )|a;< 0 }, = {( 0 , |/) 1 > 0 }, E^ = R^^\Je\ 

i=l 

The event sets are neither open nor closed. We dehne associated vector 
helds P, j G 2 *, on by 

fix, y) = (1, 0), fix, y) = (0, -1), f = f + f. 

It is a simple exercise to verify that the network is amenable and proper 
but that the associated semiflow $ : x M+ —)■ is not continuous 
along E^ or E"^ (it is continuous at ( 0 , 0 )). 

(2) Suppose that the event set E^ is the cusp dehned by {ix,y) G 

I a; 7 ^ 0, y‘^ = and \ E^. In this case any smooth 

(C^ suffices) vector held on which is tangent to E^ must vanish at 
{(0,0)} (an example of such a vector held is fax. Say), a G M). If 
we require amenability, then all traiectories which meet E^ will never 
leave E\ 0 

Proposition 4.25. An amenable asynchronous network is proper. 

Proof. We give details for the case when M is compact. Fix X G M. 
Suppose that (fi • [0, sf) —)■ M are forward trajectories for F through 
X, i G 2. Using uniqueness of solutions of diherential equations and 
dehnition 14.221 21. it is easy to see that = 02 on [0,si) fl [ 0 ,S 2 ). It 
follows that if we dehne 

T = sup {t I there is a trajectory f : [0, f) —?■ M through X} 

then we have a unique trajectory 0 : [0,T) —)■ M through X. If 
T = 00, we are done. But if T < 00, then we can extend 0 to [0,T] 
by 0(T) = lim^^T-0(0 fremarks 14.231( 3)1. If 0(T) G E°‘ then by 
dehnition I4.22lf 2). 0 extends to [0,T -|-f(0(T))), where f(0(T)) > 0. 
This contradicts the maximality of T and so T = 00 . □ 

Remarks 4.26. (1) Proposition 14.251 says nothing about the number 
of changes in the event map that occur along a trajectory. Without 
further conditions, there may be a countable inhnity of changes with 
countably many accumulation points (see dehnition 14.151 and note the 
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analogy with Zeno-like behavionr m)- 

(2) As shown in examples I4.24lf l). the semiflow given by proposi¬ 
tion |1]25] need not be continnous (as a fnnction of (X, t)). 

(3) Amenability is snfficient bnt not necessary for properness. 

4.8. Semifiows for amenable asynchronous networks. Assnme 
94 is an amenable asynchronons network with network vector held F. 
For each a & A, denote the flow of f“ by <F". 

Let X G M and (p : M+ —)■ M be the maximal integral cnrve throngh 
X for F. If follows from the definition of integral cnrve and amenability 
that there is a conntable closed snbset D = D(X.) of M+ U {cxo} snch 
that for each u E D, there exist nniqne a E A, v = v{u) E D snch that 

(m, n) n = 0, £{y) ^ a, 0 ([m, n)) C 

(For £{u) = a we need amenability.) 

Proposition 4.27. Let^ be an amenable asynchronous network. Sup¬ 
pose that for all X G M, D(X.) is finite and set D{X.) = {tj^ | 0 = < 

tf- < ... < < t]^_^_l = oo}, of = £{4){tf )), j E N*. The semiflow 

$ : M X R_|_ -E M for F is given in terms of the flows by 

$x(t) = <!>“?(••• <F“"(<I><(X, tf),tf-tf)---,t- tf), 

where t E [tf,tf+i), p E N*. 

Proof. For t E (Xp,f) is the solntion to X'(f) = f“? (X) 

with initial condition Xp = <Fx(t^). □ 

4.9. Asynchronous networks with additive input structure. A 

natural source of asynchronous networks comes from synchronous net¬ 
works with additive input structure. The event map can be either state 
dependent (with constraints) or stochastic (see the following section). 

Fix a fc node synchronous network J\f with additive input structure 
and network vector held f = (/i,..., /fc) given by. 

ei 

(4) /i(xi; Xj,,..., Xj^J = Fi(xi) -f Fy,(xj,, Xi), i E k. 

s=l 

On account of the additive input structure, it is natural to remove and 
later reinsert connections between nodes. 

For i G k, let (IFi, Cfl be the constraint defined by the 0-dimensional 
foliation of Wi = Mi. If dynamics on Ni is constrained, then dynamics 
is stopped: x' = 0. Let F be the network graph determined by (jl]) 
with associated 0-1 matrix 7 G M{k). Take P = (!,...,!) and let 
A C M,{k) be a generalized connection structure such that 
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(1) (0|7) e A 

(2) for all a = | a*') the matrix defines a subgraph of F, and 

(3) Q!jo £ {0,1} for alH G k, a G A. 

For each a & A, dehne the a-admissible vector held f“ by 


/r(xAX 


■311 


1 ^ 3 . 


.) = (!- Oio) F;(xi) + f G k, 


s=l 


and set = {f“ | a G A}. If we choose an event map £ ■. yi ^ A 
and take T = {f“ | a G A}, then 91 = (A/", A,J^,£) is an asynchronous 
network. We refer to 9t as an asynchronous network with additive input 
structure. 

For a G Al, i G k, let J{i,a) = {j \ aij = l,j G k*} be the depen¬ 
dency set of /". 


Definition 4.28. An asynchronous network 91 is input consistent if for 
any node W and a,(d & A with dependency sets satisfying J{i,a) = 
wehave/“ = /f. 

As an immediate consequence of our constructions we have 


Lemma 4.29. Asynchronous networks with additive input structure 
are input consistent. 


In summary, if 91 is an asynchronous network with additive input 
structure all the admissible vector helds are derived from the network 
vector held of a synchronous network. 


4.10. Local clocks on an asynchronous network. In this section 
we describe local clocks on an asynchronous network. We give only 
brief details sufficient for the examples we give later (the general set 
up appears in [15 )• Roughly speaking, a local clock will be associated 
to a set of nodes, or connections, and may be thought of thought of 
as a stopwatch with time r G M+. In particular, the local clock will 
run intermittently and switching between on and oh states will be 
determined by thresholds. 

Fix a hnite set of nodes A/" = {A^"o, A^i,..., N^} with associated phase 
spaces Mi, z G k, a generalized connection structure A C M,{k) and a 
constraint structure C. Local clocks will be dehned in terms of strongly 
connected components of elements of A. 

Suppose that a E A and let /9,7 be distinct strongly connected 
components of a with respective node sets A C k, R C k*. A local 
time G M+ will be dehned on f3 (or the nodes A) if there exists a 
connection Nj ^ Ni, j E B, i E A. 
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Examples 4.30. (1) The constraining node Nq is always a strongly 
connected component of a. If a = iVo —)■ iVj, then we may take {3 = 
7 = {iVo} and dehne the local time on iVj. 

(2) If a = iVo —)■ iVj -H- Nj Nq, then we may take (3 = Ni Nj, 
7 = {iVo} and obtain the local time rg = defined on Ni,Nj (or 
Ni 77 Nj). 0 

Choose a set ti, ... ,Ts of local times and set 

7" = = {r = (n,..., Ts) I Ti,..., Ts G M+}. 

We extend the phase space of A/” to Ad = M x T. Given a G A, an 
a-admissible vector held f" on M will be a smooth vector held of the 
form 

r(X, r) = (/f (X, r),..., /,“(X, T),h,...,hs), 

where hi,... ,hs ^ {0,1} are constant vector helds. 

Just as before, we dehne an Al-structure J^, an event map T ; AJ —)■ 
A and associated asynchronous network {N',A,3F,£). Our previous 
dehnitions and results continue to apply. 

Example 4.31. Suppose k = 1, M = {Nq, iVi}, and Mi = M. Choose 
a smooth vector held / : R ^ M such that 1 > f{x) > 0 for all 
a; G R. Dehne A = {0, a = Nq ^ Ni}. Dehne the local time r G R+ 
associated to a. Set Ad = R x R+. Dehne A = {f®, f"} by 

f0(a;, r) = {f{x), 0), f"(a;, r) = (0,1), {x, t) e M. 

Fix T > 0 and dehne the event map T : Ad —)■ Al by 

S{x,t) = 0, if a; 7 ^ 0 or r > T 
= a, if a; = 0 and r <T 

The asynchronous network {N',A,3F,8) is amenable. If we initialize 
at (a;o,0), xq < 0, then the system evolves until a; = 0, stops for local 
time T seconds and then restarts. In practice, the local clock is reset 
to zero after the system restarts. <0 

4.11. Stochastic event processes and asynchronous networks. 

Given node set A/”, constraint structure C, generalized connection struc¬ 
ture Al and Al-structure 3F, an event proeess is a state dependent sto¬ 
chastic process £(t,X} faking values in Al. 

Definition 4.32. (Notation as above.) A stochastie asynchronous net¬ 
work 91 is a quadruple {N',A,3F,£), where £ = £{tx) is an event pro¬ 
cess. 
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In the most general case there are no restrictions on the process 
£(t,-x.y there may be (stochastic) dependence on time t G M’*', pure 
space dependence (i^p,x) = or both. If is independent of 

time, then the event process reduces to an event map £ \ Is/l ^ A. If 
£{tX) is independent of X, then under mild conditions on £, such as 
assuming £ is Poisson, integral curves on the stochastic asynchronous 
network (A/", Al, -F, £t) will be almost surely piecewise smooth. 

We discuss stochastic asynchronous networks in more detail in [T3]. 
We give one simple example here related to additive input structure. 

Example 4.33. We follow the assumptions and notational conventions 
of section 14.91 and assume given a synchronous network with additive 
input structure and dynamics given by (|1]). Let A he & generalized 
connection structure and £ he a time dependent event process taking 
values in A. Assume M is compact and the set of times fo < < • • • 

where the connection structure changes has Poisson statistics. The 
stochastic asynchronous network {Af,A,J^,£) is an example of a sto¬ 
chastic asynchronous networks with additive input structure. Almost 
surely, trajectories will be piecewise smooth and dehned for all positive 
time. <0 


5. Model examples of asynchronous networks 

In this section, we describe two asynchronous networks using the 
formalism and ideas developed in the previous section. We refer also 
to [15 , for the detailed description of an asynchronous network mod¬ 
elling spiking neurons, adaptivity and learning (STDP). 

5.1. A transport example: train dynamics. We use a simple trans¬ 
port example - a single track line with a passing loop - to illustrate 
characteristic features of asynchronous networks in a setting requiring 
minimal structure and background knowledge. 

Consider two trains Ti,T 2 travelling in opposite directions along a 
single track railway line; see hgure [3l We assume no central control 
and no communication between train drivers unless both trains are in 
the passing loop. 

Take as phase spaces for the trains the closed interval I = [—a,b], 
where a, 6 > 0. Suppose the end points of I correspond to the stations 
A (at —a) and B (at b) and that the passing loop is at 0 G I. Assume 
that the passing loop is associated with a third station P. 

The position of train Tj at time t > 0 will be denoted by Xi(t) G /, 
i G 2. Suppose that a:i(0) = —a, X2(0) = b. Assume that, outside 
of the stations A, B, P, the velocity of the trains is given by smooth 



ASYNCHRONOUS NETWORKS 


29 


A P B 



“II ‘12 


Figure 3. Two trains on a single track railway line with 
a passing loop and stations. 


vector fields Fi, V 2 : / —)■ M satisfying 

Vi{x) > 0 > V 2 {x), X e /. 

That is, Ti is moving to the right and T 2 to the left. In order to pass 
each other, the trains must enter the passing loop and stop at P. 

Fix thresholds S, S'!, S' 2 , Ti, T 2 G M+. Train % will depart at time Ti, 
i G 2. We require that trains have to be together in station P for time 
S and, additionally, the train Tj must be in the station for time Si, 
i G 2 (this is an additional condition on T* only if Si > S). The trains 
can move out of the station when these thresholds are met. Note that 
the trains will not generally leave the station at the same time if Si > S 
or S '2 > S. We model train dynamics by an asynchronous network. 

First we discuss connection structures. Associate the node W with 
train Tj, i G 2. Train Tj will be stopped at P only if there is a 
connection «* = A"o —)■ Ni, i G 2. We only allow communication 
between trains when both trains are stopped at P. In this case, the 
connection structure will be /3 = A"o —)■ W t-G A "2 t— A" 2 - If either train 
is not stopped at P, there is no connection between the trains. 

As the drivers of the trains cannot communicate (unless both trains 
are in the station P) and there is no central control, the times associ¬ 
ated with the thresholds Si, S 2 will be local times. Specifically, when 
train %i stops at P, the driver’s stopwatch will be started. This will 
be a local time r* for Tj and associated to the connection Nq —)• W, 
When both trains are stopped at P, we use a third local time r = ri 2 
associated to the connection Ni gg N 2 (alternatively, the drivers could 
synchronize their stopwatches but still the stopwatches may not run at 
the same speed). 

We describe this setup using our formalism for asynchronous net¬ 
works. As network phase space we take 


M = {(X,r) = (xi,X2,ri,r2,r) | Xi,X2 G I, Ti,T2,t G M+} = /^ x 
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We define the generalized connection structure A = {ai, 0 : 2 , (3, 0} and 
let T be the ^-structure given by 

f®(X,T) = ((W(Tl),l^2(a;2)), (0,0,0)) 
rHX,T) = ((0,l^2(a;2)), (1,0,0)) 

P^(X,t) = ((W(a;i),0), (0,1,0)) 
f^(X,T) = ((0,0), (1,1,1)) 

We dehne the event map £ ■. M. Ahy 

'ai if (xi = 0, X2 > 0) V ((xi = 0, X2 < 0) A (ri < 5i)) 

^ ^ 012 if (x2 = 0, xi < 0) V ((x2 = 0, xi > 0) A {t2 < S 2 )) 

13 if (xi = X 2 = 0) A ((r < 5) V ((ri < Si) A {t 2 < S 2 ))) 

0 otherwise. 

Here we have used the logical connectives V for or and A for and. 
Dynamics on the asynchronous network 91 = {N',A,3F,8) is given by 
the vector held F(X) = f^*^^^(X). Provided that we initialize so that 
a;i(0) < 0 < 0 : 2 ( 0 ), ri(0) = T2(0) = r(0) = 0, it is easy to see that 91 is 
amenable. 

5.1.1. Initialization, termination and function. The network 91 has a 
function: each train has to traverse the line to reach the opposite sta¬ 
tion. Thus we can regard 91 as a functional asynchronous network. 
Formally, dehne initialization and termination sets by Ii = {—a}, 
I 2 = {h} and Fi = {6}, F 2 = {—«} respectively. We call I = Ii x I 2 and 
F = Fi X F 2 the initialization and termination sets for 91. The function 
of the network is to get from I to F in hnite time. 

Typically, the thresholds S', S'!, 5*2, Ti, T 2 G M+ will be chosen stochas¬ 
tically. For example, the starting times Ti,T 2 according to an ex¬ 
ponential distribution. If we initialize at (—a, Ti), (6, T 2 ), and take 
ri(0) = T 2 ( 0 ) = r(0) = 0, it is easy to verify that solutions will be 
dehned and continuous for all positive time under the assumption that 
a train stops when it reaches its termination set. 

5.1.2. Adding dynamics. The trains only “interact” when both are 
stopped at P. We now add a non-trivial dynamic interaction when 
the trains are stopped at P. To this end, we additionally require that 

(1) The drivers are running oscillators of approximately the same 
frequency (randomly initialized at the start of the trip). 

(2) When both trains are at P, the oscillators are cross-coupled 
allowing for eventual approximate frequency synchronization. 

(3) The trains cannot restart until the oscillators have phase syn¬ 
chronized to within e, where 0 < e < 0.5. 
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For example, fix M and define Hiff) = A;sin2710, 0 G T, where 

/c > 0. Take as network phase space x T^. Dehne vector 

helds h® = and on Al* by 

h®(X,r,0i,02) = (0,0,a;i,a;2) 

h^(X, T, 01, 02) = (0, 0,071 + i/(02 — 0i), 072 + -fA(0i — 02)) 

Dehne a new A-structure A* by 

g® = f® + h®, = pi + h“b g“2 = + h"^ g^ = + h^, 

where G A do not depend on (0i,02) G T^. Modify the 

event map E by requiring that £^(X, r, 0i, 02) = /5 iff 

{xi = X 2 = 0) A ((r < S') V (|0i - 02| > e) V ((ri < Si) A (r 2 < S' 2 ))) 

In this case, for almost all initializations, the oscillators will eventually 
phase synchronize to within £ provided that sin“^(|o7i — 072 I/ 2 A;) < 
27ie. In particular, if 07i = 072 , the oscillators will synchronize unless 
|0i(O)-02(0)1 = 0.5. 

5.1.3. Relations with Filippov systems. Assume all the thresholds of 
our model are zero. Note that if S' = S'! = S '2 = 0, then there is no need 
for local clocks and we may model by the asynchronous network DT* = 
(AT, A*, A*,S*), where A* = {q!i,Q! 2,0}, A* = {f®,Pbf “^}5 where 
f®(X) = {Vi{xi),V2{x2)), P^X) = {0,V2{X2)), P^(X) = (Pi(a:i),0), 
and the event map E* is dehned by 

{ «! if Xi = 0, X 2 > 0 
02 if a;2 = 0, xi < 0 
0 otherwise. 

We show dynamics for 91* in hgure H] under the initialization assump¬ 
tion that Xi(0) < 0 < ^ 2 ( 0 ). Referring to the hgure, trajectory rj 
corresponds to train T 2 reaching P hrst and restarting only when Ti 
reaches P. Train Ti reaches P hrst for the trajectory u. Regardless of 
which train reaches P hrst, the ‘exit trajectory’ cj) is always the same 
and so there is a reduction to 1-dimensional dynamics. If both trains 
arrive simultaneously at P, neither stops. 

The dynamics shown in hgure 0] is suggestive of a Filippov system [271 
E] and it is natural to ask whether there are connections between asyn¬ 
chronous network and Filippov systems. Set = {{xi,X 2 ) \ X 1 X 2 < 0} 
and observe that dynamics on 91* is given by a continuous semihow 
$* : Rg X R+ —>■ Rq. We dehne a Filippov system on R^, with contin¬ 
uous semihow $ : R^ x R+ —>■ R^, such that d> = <F* on R^. To this 
end we let Qp-, i,j G {+, —} denote the closed quadrants of R^ (so 
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Figure 4. Dynamics on a one track line with passing loop. 


*2 



Figure 5. Dynamics for the Filippov system. Trajec¬ 
tories 1 ] and (j) are nnchanged; trajectories k and ^ corre¬ 
spond to one train reversing after the other train enters 
the passing loop and are artifacts of the Filippov repre¬ 
sentation. 


Q+_ = {{xi,X2) I Xi > 0,X2 
on each qnadrant by 

V++(a:i,a;2) = 
V+_(a;i,a;2) = 
V__(a;i,a;2) = 
Y_+{Xi,X2) = 


< 0}, etc) and define smooth vector fields 

i-V{Xi),V2{x2)), iXi,X2) e Q++ 
iV{Xi),V2{x2)), {Xi,X2) e Q+- 
{V{Xi),-V2{x2)), {Xi,X2) e Q— 
(y{Xi),V2{x2)), {Xi,X2) e Q-+. 


These vector fields nniqnely define a smooth vector field V on the 
nnion of the interiors of the qnadrants. We extend V to a piecewise 
smooth vector field on \ {(0)0)} using the Filippov conventions. 
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Thus, we regard the Xj-axis as a sliding line S\ i G 2, and define V 

on n dQ _= 5*1 to be the unique convex combination of 

V__|_ and V_which is tangent to (in this case (V_+ V_)/2). 

Finally define V(0,0) = (Vi(0), 1 / 2 ( 0 )). The piecewise smooth vector 
field V has a continuous flow <F : x R_|_ —)• (integral curves are 

defined using the standard conventions of piecewise smooth dynamics 
- see [27]) and <F|Ro = d>*. Of course, the semifiow on R^ \ R^ does 
not have an interpretation in terms of trains on a line with a passing 
loop (see figure Ej) • 

In an asynchronous network, dynamics on event sets is given explic¬ 
itly rather than by the conventions used in Filippov systems. However, 
as we have shown, asynchronous networks can sometimes be locally 
represented by a Filippov system (see [T2| for more details and greater 
generality). This relationship suggests the possibility of applying meth¬ 
ods and results from the extensive bifurcation theory of nonsmooth 
systems to asynchronous networks. 

5.1.4. Combining and splitting nodes. We conclude our discussion of 
asynchronous networks modelling transport with a brief description of 
processes defined by combining or splitting nodes (a dynamical version 
of a Petri Net na). We consider the simplest cases of two trains com¬ 
bining to form a single train or one train splitting to form two trains. 
We only give details for the first case but note that both situations are 
easily generalized and also, like much of what we have discussed above, 
apply naturally to production networks. 

Consider node sets J\f°‘ = {Nq, Ni, N 2 } and = {iVo,iVi 2 }, where 
Ni, N 2 , Ni 2 have phase space R and correspond to trains ‘Xi,T 2 , T 12 
respectively. We give a network formulation of the event where trains 
Ti, 22 are combined to form a single train T 12 (see figure (Hj). Fix vector 


Direction of time and space evoiution 



Figure 6. Combining two trains into a single train. 
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fields Vi,V 2 ,Vi 2 on M and assume Vi(x), V 2 (x), Vi 2 (a;) > 0 all x G M. 
Define generalized connections structures 

= { 0 , Oil = Nq Ni, a2 = Nq —)■ N 2 , = Nq Ni -H- N 2 ^ -^ 2 }, 

= {0,7 = iVo^iVi 2 }. 

Assume a local clock with time r = ri 2 that is shared between the 
connection (3 G and 7 G A^. Dehne network phase spaces for 
J\f^ to be X M+, Ai^ = M x M+ respectively. Define the 

^“-structure by 

r .0 = 14 ), 0 ), = (( 0 , D 2 ), 0 ), = (lA, o), = ((o, o), i). 

and the ^^-structure by = (V 12 , 0 ), = ( 0 , 1 ). 

Fix thresholds S 2 , Si > 0. The threshold Si gives the time taken to 
combine Ti and T 2 , and S 2 models the time T 12 spends in the station 
before leaving. Initialize J\f°‘ so that Xi( 0 ),X 2 ( 0 ) < 0 and r(0) = 0. 
The event map r) is defined for Xi, X 2 <0 and r < S' by 

S“(X,r) = 0, Xi,X 2 < 0 

= Di, Xi = 0, X 2 < 0 
= 02, Xi < 0, X 2 = 0 

= /3, Xi = X 2 = 0,T < Si 

The event map S^(xi 2 , r) is dehned for X 12 > 0 and r > Si by 

S'’(xi 2 ,r) = 7 , X 12 = 0 ,r < Si + S 2 
= 0, otherwise 

When T = Si, we switch from network J\f°‘ to Af’. 

The splitting construction is similar except that we need to split 
the local clock for the combined train into two clocks, one for each 
separated train. 

5.2. Power grids and microgrids. 

5.2.1. Power grids as asynchronous networks. We hrst consider an un¬ 
realistic, but simple and instructive model that shows how asynchro¬ 
nous and event dependent effects can naturally £t into the framework 
of power grids. In the following section, we describe how more realis¬ 
tic models are obtained, their limitations, and where we might expect 
asynchronous network models to be useful. 

We use the simplest model [28] for power grid frequency stability that 
assumes generators are synchronous, loads are synchronous motors and 
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consider the network of mechanical phase oscillators 

n 

(5) e'- + ajO'j = hj sin( 6 'i - Oj), j G n, 

i=l 

where {kij) is a a symmetric matrix, all entries positive (zero is al¬ 
lowed). If Ylj Pj ~ 0 system can reach an eqnilibrinm {Pj < 0 
corresponds to a load). Let T be the (nndirected) graph determined by 
the matrix of connections given by {kij). While the network described 
by (Ej) is not asynchronous (and the main interest lies with the stabil¬ 
ity of the equilibrium solution), the dynamics of real-world power grids 
are subject to factors that cannot be adequately described by a syn¬ 
chronous model. For integrity of transmission lines, as well as system 
stability, it is essential that the phase differences 16^4 — 6j\ are bounded 
away from 7 r/ 2 . For example, we might require \6i — 6j\ < Tij, where 
Tij G (0, vr/2) will be a threshold determining the safe operational load 
for the transmission line. This leads to the construction of state depen¬ 
dent event maps £ij : — )■ {F,F \ {z GG j}}. If \9i — 9j\ > Tij, then 
s,,{e) = F \ {z GG j} and the transmission line between nodes z and j 
is disconnected. Equation ([5]) is modihed accordingly. Similarly, lines 
or generators may be disconnected because of external events - such as 
lightening strikes or mechanical breakdowns. These can be modelled 
using a stochastic event map. 

As indicated above, this model is unrealistic (it is not true, for ex¬ 
ample, that typical loads are synchronous motors). In the next section, 
we indicate how more realistic models are obtained, their limitations, 
and where we might expect asychronous network models to be useful. 

5.2.2. Network-reduced model for power grids. We give an overview of 
the network-reduced coupled phase oscillator model for power grids, 
largely based on Dorfler [20] , and refer the reader to [53l |20| for greater 
generality, alternative models, and the many details we omit. Apart 
from describing the model, our goal is part cautionary (it is not evident 
that general theories of synchronous or asynchronous networks have 
much to contribute to stability problems involving structural change), 
and part comparative with the models we describe later for microgrids. 

Assume a power grid with synchronous generators, DC power sources, 
transmission lines and various types of load. We assume a reference 
frequency ojr for the power grid, usually 50Hz or 60Hz, and note that 
frequency synchronization is critical for the stability of the power grid: 
our equations will be written nominally in terms of phases 9i{t) but 
for the models, we can always replace 9i{t) by 9i{t) — u^t to get the 
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(same) equations for phase deviations that are needed for stability the¬ 
ory (phase differences, but not absolute phases, matter). 

Formally, assume given an undirected (connected) weighted graph Q 
with node set V = n and edge set S gV^. Nodes will be partitioned as 
V = Vi U V 2 U V 3 , where Vi consists of synchronous generators, V 2 are 
DC power sources, and V 3 comprises various types of load (see below 
and note we do not consider all types of load). 

Each edge {i,j) G £, i 7 ^ j, is weighted by a non-zero admittance 
Yij G C and corresponds to a transmission line. The imaginary part 
^(Yij) is the susceptance of transmission line and iH(hij) is the conduc¬ 
tance. Typically, a high voltage AC transmission line is regarded as 
lossless {yi{Yij) = 0) and inductive {'3{Yij) > 0). We allow self-loops 
i = j, these will correspond to loads modelled as impedances to ground 
(nonzero “shunt admittances”). 

To each node is associated a voltage phasor V) = correspond¬ 

ing to phase 9i and magnitude \Vi\ of the sinusoidal solution to the 
circuit equations. 

For a lossless network, the power flow from node i to node j is given 
by ajj sin( 6 'j — 0 ^), where aij = |V)| |f^|J(y)j) gives the maximal power 
flow (see Kundur [101 Chapter 6 ]). 

5.2.3. Synchronous generators. We assume dynamics of synchronous 
generators are given by 

n 

(6) Mi9" -\- Di9[ = Pm,i + up sin( 0 j — 9i), i E Vi, 

i=i 

where 9i,9^ are generator rotor angle and frequency, Mi,Di > 0 are 
inertia and damping coefficients, and Pm,i is mechanical power input. 

5.2.4. DC/AC inverters: droop controllers. Each DC source in V 2 is 
connected to the AC grid via a DC/AC inverter following a frequency 
droop control law which obeys the dynamics [ 00 ] 

n 

(7) Di9[ = Pd,i + ^ Oij sm{9j - 9i), ie V 2 . 

i=i 

5.2.5. Frequency dependent loads. We assume the active power demand 
drawn by load i consists of a constant term Pi^i > 0 and a frequency 
dependent term Di9[, Di > 0, leading to the power balance equation 

n 

Di9[ = -Pi^i + ^ aij sm{9j - 9i), i E V 3 ,/, 
j=i 


(8) 
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where V 3 j is the subset of V 3 consisting of frequency dependent loads. 
Equation ([H]) is of the same form as ([7]), and we may replace V 2 by 
V 2 U V 3 j and consider the general equation 

n 

( 9 ) Di9[ = Wj + ^ Qij sin( 0 j -9i), i e V 2 , 

i=i 

where Ui is positive if the node is a DC generator and negative if it is 
a frequency dependent load. 

We can similarly allow for loads which are synchronous motors, in¬ 
corporate them in Vi and consider 

n 

( 10 ) Mi9” -\- Di9[ = L/Ji -\- ^ ^ Qij sin[9j — 9i), i G Vi, 

i=i 

where cuj is positive if the node is a synchronous generator and negative 
if it is a synchronous motor. 

5.2.6. Constant current and constant admittance loads. We assume the 
remaining loads each require a constant amount of current and have a 
shunt admittance (to ground). In this case we have a current balance 
equation and, through the process of Kron reduction [22] , may obtain 
a reduced network the equations of which are 

n 

( 11 ) Mi9” -\- Di9[ = u:i + ^ ^ dij sin[9j — 9i (fij), i G Vi, 

i=i 

n 

( 12 ) Di9'- = Dj -f dij sin(6'j — 9i + ^ij), i G V 2 . 

j=i 

We refer to [2l| for the explicit form of the coefficients in (11111121) . 

The original power grid network is typically sparse with many nodes 
- V 3 is large. The process of Kron reduction results in a much smaller 
network which will be all-to-all coupled provided that the graph dehned 
by V 3 is connected [22]. However, even if the original transmission 
lines are lossless, the phase shifts (pij will generally be non-zero and not 
necessarily always small (we refer to [521 §6-2 Figure 4] for data from 
a real power grid network). The presence of phase shifts can and does 
make it harder to frequency synchronize ( 11111121 ) . 

From the point of view of transmission line failure in a power grid, 
even if the removal of an edge still results in a all-to-all coupled re¬ 
duced network, many of the coupling coefficients dij will change. It is 
a hard problem, that goes beyond existing analytical theory for syn¬ 
chronous and asynchronous networks, to get good insight into whether 
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or not a breakdown will destabilize the network (this is irrespective of 
phenomena like Braess’s paradox [Oil El])- 

5.2.7. Microgrids. Assume given a stable power grid network, robust 
to “small” changes in power demand, and consider the problem of 
modelling a microgrid and its combination or separation from the main 
grid. We outline structural and logical issues to make transparent the 
connection with asynchronous networks and largely ignore dynamics so 
as to keep the model simple and our discussion short (we refer to [231 
mMM\ for more details and references on microgrids and control 
from a large and rapidly developing literature in this area). Assume 
power generation in the microgrid is from DC generators (such as solar 
power or DC wind power) and that Vi = 0 (most motor loads are not 
synchronous). Assume the microgrid is Kron reduced. 

Unlike the power grid model described above, we allow directed (one 
way) connections and a constraining node. Consider the simplihed net¬ 
work Af = {Nq, Nb, Ng, Np}, where the nodes Np, Nq, Np correspond 
to a large capacity battery (buffer), a DC generator, and main power 
grid respectively, and dehne subnetworks A/m = {Np, Nq} (microgrid) 
and A/p = {Np} (main power grid). 

The battery acts as reserve storage or buffer for the microgrid; in par¬ 
ticular to maintain power in the event of intermittent loss of generated 
DC power or when the microgrid has been separated “islanded” from 
the main power grid. We suppose battery capacity B = B(t) E [0, Bm], 
where Bm corresponds to the battery being fully charged. We suppose 
that the DC generator produces power O = 0(t) G [0, Om], where Om 
is the maximum power than can be generated. 

The constraining node will play a role when the microgrid is islanded 
and is to be reconnected to the main power grid: either because the 
microgrid has insufficient power for the microgrid loads or because the 
microgrid has an excess of available power some of which can now be 
contributed to the main power grid. In either case a transition process 
needs to be implemented where the droop controller for the DC/AC 
converter needs to bring the AC output of the microgrid in precise 
voltage (phase, frequency and magnitude) synchronization with the 
state of the power grid at the connection point (s) to the microgrid. 
Similarly, we can constrain when the microgrid is to be islanded from 
the main grid so that the reduction in power contributed to the main 
power grid is gradual and done over an appropriate time scale so as 
not to destabilize the main power grid. 

Leaving aside the dynamics of islanding and combining the microgrid 
with the main power grid, the generalized connection structures and 
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control logic we need for management of the microgrid are complex and 
depend on several thresholds which may need to be time dependent 
- for example, if we use a time dependent model for the projected 
microgrid power load. If the microgrid is islanded, we work with Mm 
and use the generalized connection structure 

AM = {a = NG^ Nb, /3 = Nb^ Ng, 0 }. 

The connection structure a corresponds to the DC generator having 
sufficient output to supply all power needed for the microgrid load and 
with a surplus which can be used to charge the battery, (3 corresponds to 
battery and generator providing all necessary power for the microgrid, 
and 0 corresponds to the generator providing all needed power for 
the microgrid and either there is surplus power available for battery 
charging or the battery is fully charged. Thresholds that determine 
switching between these states are chosen so as to avoid “chattering” 
in the control system. 

If the microgrid is combined with the main power grid, this can be ei¬ 
ther because battery and DC generators cannot provide sufficient power 
for the microgrid load or because the microgrid has surplus power which 
can be contributed to the main power grid or because the main power 
grid is stressed (possibly locally detected by frequency variation) and 
the battery state of the microgrid is sufficiently high to allow a tempo¬ 
rary power contribution to the main grid. As generalized connection 
structure A we take the set of connection structures 

Ng Nm, Ng —)■ Nm t— Nb, 

Nb Ng Nm, Nm —t Ng, Nm —t Ng —)■ Nb, 

Each of these connection structures has a natural interpretation. For 
example, Nm —t Ng —)■ Nb corresponds to the main power grid con¬ 
tributing to both the load of the microgrid and battery charging while 
Ng —)■ Nm t— Nb means battery and DC generator are contributing 
power to the main power grid as well as supplying all the power for the 
microgrid. On the other hand, Ng —)■ Nm means DC generated power, 
but not battery power, is being contributed to the main power grid. 

Of course, what we have described above is highly simplihed as we 
have taken no account of (1) multiple DC generators and batteries 
within a microgrid, or (2) multiple microgrids. In the latter case, we 
need to take care that microgrid switching does not synchronize as this 
could lead to large destabilizing changes in load on the main grid. 
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6. Products of asynchronous networks 


We conclude part I with the definition of the product of asynchronous 
networks and give sufficient conditions for an asynchronous network to 
decompose as a product of two or more asynchronous networks. Al¬ 
though the methods we use are elementary, the study of products is 
illuminating as it clarifies some subtleties in both the event map and 
the functional structure that are not present in the theory of synchro¬ 
nous networks. These ideas play a central role in the proof of the 
modularization of dynamics theorem in part II. 


6.1. Products. Given a,/3 G M{k), define a\J (3 G M{k) (the join 
of a and (3) by 

(a V (3)ij = ma.yi{aij,(3ij] , b J ^ k 

(the max-plus addition of tropical algebra |35])- have a V 0 = a for 
all a G M{k). \i A, B C M{k) are generalized connection structures, 
define the generalized connection structure Ay B hy 


Ay B = {a\J I3 \ a & A, (3 E B} . 

Note that 0 G AyB if and only if 0 G AhB. Consequently, if 0 G AyB, 
then A,B C Ay B. 

Suppose that A is a nonempty subset of k containing kA elements. 
There is a natural embedding of M{kA) in M{k) defined by relabelling 
the matrices in M{kA) according to A. Specifically, map the matrix 
{aij)ij^A £ M{kA) to the matrix a E M{k) defined by 


aij — 



for z, j G A, 
otherwise. 


This embedding extends to an embedding of M,{kA) in M,{k) by 


DiO — 



for i E A, 
otherwise. 


Given disjoint nonempty subsets A, i? of fc, regard M,{kA), M,{kB) as 
embedded in M,{k). Given a E M,{kA), [3 E M,{kB), define 

ay (3 = a y /3 E M,{k). 

This extends to the join Ay B of generalized connection structures on 
disjoint sets of nontrivial nodes. 

Let Af = {A"o,..., Nk} and A be a proper subset of k. Define = 
{Nj\j G A*} and = HiGT Denote points in Myi by 
Suppose P = k \ A. We have = {Xq} and x M. 

If C'^, O'® are constraint structures on A/"^, respectively, let C = 
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V denote the induced constraint structure on Af - well defined 
since constraints depend only on nodes and An B = ^. 

More generally, given disjoint node sets = {Nj \j E A*}, Af^ = 
{Nj \j E B*}, we can identify A^B with complementary subsets of k, 
where k is the total number of elements in AU B^ and then follow the 
conventions described above. 

Definition 6.1. (Notation and assumptions as above.) Given asyn¬ 
chronous networks X E {A,!?}, define the 

product X to be the asynchronous network 91 = {Af,A,X,S) 
where 

( 1 ) Af = Af^uAf^, 

(2) C = V C^, 

(3) A = A^yA^, 

(4) X = X E A^, E A^}, and 

(5) S is defined by 

S{Xa, Xb) = S^{Xa) V S^{Xb), for (X^, X^) eMaX M^. 

Remark 6.2. If 91"^, 91^ are proper (or amenable), then so is 91"^ x 91^. 

Lemma 6.3. (Notation of definition \6.1\ ) The network vector field on 
91"^ X 91'® is given by 

(13) F(X^,Xb) = (ff (''")(X^),ff (^"^Xs)), 

for all (XaiXb) E IVfyi x IVI^. 

Proof. Immediate from the definitions. □ 

6.2. Decomposability. 

Definition 6.4. An asynchronous network {Af,A,fF,£) is decompos¬ 
able if it can be written as a product of asynchronous networks. If the 
network is not decomposable, it is indecomposable. 

Example 6.5. Suppose that TV is a synchronous network with con¬ 
nection structure a E M{k) and o-admissible network vector field f 
satisfying conditions (Nl-3) of section [2l Since a encodes the depen¬ 
dencies of f it is trivial that Af can be written as a product of two 
synchronous networks iff the network graph r„ is disconnected. <0 

Our aim to find sufficient conditions on an asynchronous network for 
it to be decomposable. 

Definition 6.6. The connection graph of the asynchronous network 
91 = {Af, A,J-,8) is the graph defined by the 0-1 matrix Tg^ = VaeA ■ 
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Lemma 6.7. If an asynchronous network is decomposable, then the 
connection graph T<yi of ^ has at least two connected components. 

Proof. If fit is decomposable, then DT = where A, B are proper 

complementary subsets of k. Since there are no connections between 
nodes in Af^ and Af^, has at least two connected components. □ 

Remark 6 . 8 . Lemma [6.71 gives a necessary condition for decomposabil- 
ity which is not sufficient. There are two issues. First, the event map 
encodes information about spatial dependence of node interactions that 
cannot be deduced from the connection graph. Second, the admissi¬ 
ble vector fields may have dependencies that are incompatible with 
decomposability. 


Example 6.9. Let k = 2, Mi = M 2 = M. Define connection struc¬ 
tures Ui = Nq ^ Ni, i E 2 and generalized connection structure 
A = {0, Q!i, 0 ^ 2 , /9 = ai V ^ 2 }. Suppose the event map is given by 


S{Xi,X2) = < 


«!, 

T2, 

/9, 

10, 


if xi < 0 , X 2 = 0 
if Xi = t),X2 > 
if Xi = X 2 = 0 , 
otherwise 


In this case, A = y Af, where A^ = {0, «*}, f G 2, and the network 
graph is disconnected. However, there is no way to write £{xi,X 2 ) as 
£^{xi) \J£‘^{x2) as the event sets involving xi E Mi depend nontrivially 
on X2 E M2. Hence the network cannot be decomposable or even 
equivalent to a decomposable network whatever choice we make for 
admissible vector fields. 

Suppose instead we define the event map by 


£{Xi,X2) = < 


«!, 

D 2 , 

/S, 

0 , 


if = 0 , a ;2 7 ^ 0 
if a ;2 = 0 , Xi 7 ^ 0 
if xi = 0:2 = 0 
otherwise 


In this case .4, = .4,^ V .4,^ and we may write £ = £^\J £‘^ where £^*(0) = 
ctj, and £'‘{xi) = 0, x, 7^ 0, i G 2. Suppose that f“i(xi,X 2 ) = ( 0 ,n 2 ), 
f" 2 (xi,X 2 ) = (ti, 0), f®(xi,X 2 ) = (ni,n 2 ), where vi,V 2 7 ^ 0. For the 
moment leave unspecified. Define IF'" = {ff,f“*}, where ff(xj) = n*, 
= 0, i G 2. Observe that f® = ff x f®, x f| and 

f “2 = f 0 X For (TV, A,IF,£) to be a product we additionally require 
f^(xi,X 2 ) = (ff^(xi), f^^(x 2 )) = (0,0), all (xi,X 2 ) G In particular, 
if f^( 0 , 0 ) 7 ^ ( 0 , 0 ), the network {Af,A,IF,£) is not even equivalent to 
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a product network. However, if f^(0,0) = (0,0), then the network 
will be equivalent to a product network if we redehne 
to be X f 2 ^ (this does not change the values of on 0 

6.3. Sufficient conditions for decomposability. Let OTbe an asyn¬ 
chronous network with k nodes and C* be a proper connected compo¬ 
nent of the connection graph Identify C with the nonempty subset 
of k corresponding to the labels of the nodes in the component C. 
Let C = k \ C*. Since C* is a connected component of Lg^, we can write 
each a ^ A uniquely as a = ac V where ac, o.-^ are connection 
structures on and respectively. Set A^ = {ac \ a G A}. We 
have a well defined projection nc ■ A ^ AP defined by 7rc(a) = ac- 
Define the event map x —)■ A^ by 

^^(Xc,X^)=7rc(^(Xc,X^)). 

Definition 6.10. An asynchronous network ffi is structurally decom¬ 
posable if for any connected component C of the connection graph 
the map is independent of X^ G (that is, £^'"(Xc',X^) = 
£^^(X( 7 ) where ; Me —)■ A^). 

Remark 6.11. Structural decomposability implies conditions on struc¬ 
tural dependencies that will generally be different from the dependen¬ 
cies of the network vector field. For example, suppose a component 
C of the connection graph contains the node Ni. If the node Ni is 
stopped there may be a condition that Ni will restart when the state 
of another node, say X 2 , attains a certain value. Necessarily, N 2 must 
he in C (structural decomposability). However, there need be no con¬ 
nection between Ni and N 2 unless C contains exactly two nodes. 

Suppose that 91 is structurally decomposable and that Fg^ has con¬ 
nected components Ci,...,Cq. Set = Mq, A^ = ttcPA), £ G 
q. By structural decomposability we may write i^(X) = 
where —)■ A^. For a G .4,, £ G q, set a^ = nePa) G A^ and 

= {8^)-\aP C M,. 

Lemma 6.12. (Notation as above.) IfOl is structurally decomposable 
and F^ has connected components Ci,... ,Cq, then 

Ea = ^ 

teq teq 

Proof. An immediate consequence of structural decomposability. □ 
If C* be a proper connected component of the connection graph F^^ 
of an asynchronous network 91, then by admissibility 

P = fg X f^, for all a G 
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where fg : Me —)• TMe and —)■ TM^. 

In order that DT be decomposable, this decomposition has to be com¬ 
patible with the projections vre : ^ t , tt^ : ^ . In particniar, 

if connections in the set of nodes that are in C are added or deleted, 
dynamics on Me is not affected. 

Definition 6.13. (Notation as above.) The asynchronons network D1 
is dynamically decomposable if for any connected component C of Tsj^, 
we have 

fOL _ f0 

'-c — '-c 

for all a^(3 & A snch that 7re(a) = Tici.fd)- 

Lemma 6.14. (Notation as above.) Input consistent asynchronous 
networks are dynamically decomposable. In particular, asynchronous 
networks with additive input structure are dynamically decomposable. 

Proof. Given i G k, a G .4, let J{i, a) be the associated dependency set 
for node Ni. a, (I E A and J{i,a) = ■I{i,(3), then /“ = /f by inpnt 
consistency. IfiEC, where G is a connected component of the network 
graph Tsj^, then J{i, a) nk C C for all a E A. Hence J(f, a) = J{i, oc V 
a-^) is independent of a-^. Inpnt consistency implies that = 

jaevj ^ ^ yields dynamical decomposability. □ 

We now state the main resnlt of this section. 

Theorem 6.15. Let ‘HI be a structurally and dynamically decomposable 
asynchronous network with connection graph T. // T has connected 
components Ci,...,Cq then there exist indecomposable asynchronous 
networks 91^, ..., such that 

m = m^ X ■■■xHlL 

Proof. For ^ G q, dehne A^ = {a£ 7r£(a) | a E A} and = {f“^ 

: M^ —)■ TM^ I a E A}. By dynamical in decomposability we have 
f" = a E A. Constraint strnctnres are dehned for 

individnal nodes and so factorise natnrally. Let : M^ —)■ A^ be the 
event maps given by strnctnral indecomposability. If we let 91^ be the 
asynchronons network where = {iVo} U {W \ i E 

Ci}, £ G q, then 91 = Heq ^ 

Onr conclnding resnlt on decomposability is an immediate conse- 
qnence of lemma 16.71 and theorem 16.151 

Corollary 6.16. A structurally and dynamically decomposable asyn¬ 
chronous network 91 is decomposable if and only if its connection graph 
has more than one nontrivial connected component. 
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6.4. Factorization of asynchronous networks. Assume for this 
section that DT = is an asynchronous network which is 

not necessarily structurally or dynamically indecomposable. 

Definition 6.17. The asynchronous network 91^ is a factor of 91 if 
there is an asynchronous network 91^ such that 91 = 91^ x 91^. 

The proof of the next lemma is immediate from the definition of a 
product. 

Lemma 6.18. //91^ is a factor o/91, then the connection graph Tg^i 
is a union of connected components ofTt^. 

Remark 6.19. If 91^ is indecomposable, the connection graph T^i may 
have more than one component - unless 91 is structurally and dynam¬ 
ically indecomposable (theorem [6T5]). 

Proposition 6.20. Every asynchronous network'^ has a factorization 
as a product of indecomposable asynchronous networks. The 
factorization is unique, up to the order of factors. 

Proof. Existence is obvious. The uniqueness of factorization follows 
easily from lemma 16.181 □ 
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